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The population of pulsars with interpulses and the 
implications for beam evolution 
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ABSTRACT 

The observed fraction of pulsars with interpulses, their period distribution and the 
observed pulse width versus pulse period correlation is shown to be inconsistent with 
a model in which the angle a between the magnetic axis and the rotation axis is 
random. This conclusion appears to be unavoidable, even when non-circular beams 
are considered. Allowing the magnetic axis to align from a random distribution at 
birth with a timescale of ~ 7 x 10^ years can, however, explain those observations 
well. The timescale derived is consistent with that obtained via independent methods. 
The probability that a pulsar beam intersects the line of sight is a function of the angle 
a and therefore beam evolution has important consequences for evolutionary models 
and for estimations of the total number of neutron stars. The validity of the standard 
formula for the spin-down rate, which is independent of a appears to be questionable. 
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1 INTRODUCTION 

Pulsars are formed from the supernova explosions which 
mark the death of high mass stars. During the pulsar birth 
process, th e pulsar's rotation axis is aligned with its veloc - 
ity vector (|Spruit fc PhinnevI ITooi : [Johnston era!] |2005| ). 
A further question arises as to the orientation of the mag- 
netic axis with respect to the rotation axis at birth and 
the evolution of this angle with time. Over the years, ob- 
servational evidence seems to support the view that the 
magneti c and rotat ion axes align over time (Candy & B laiij 
[l986 ; Lvne fc Manchester 1988 : Tauris fc Manchester 1991). 



although there have bee n claims for no evolution ( McKinnonI 
1 19931 : iGil fc Hani 119961) and even th eoretical reasons for 
counter-alignment ( Beskin et al.lll988l ). Understanding this 
issue is important, not only for theoretical considerations of 
how braking torque works in neutron stars, but also for the 
evolution of the pulsar beam. This has important implica- 
tions for evolutionary models and for pulsar surveys gener- 
ally, such as those p lanned for the Square Kilometre Array 
l|Kramer et al.ll2004l ). 

Radio emission from pulsars arises from a few hundred 
kilometres above the polar cap, the region bounded by the 
last open magnetic field lines. It can be shown simply that 
the radius of the polar cap, r, is related to the pulsar spin 
period, P, via r oc P""'^. In turn then, the opening angle of 
the beam, p, can be expressed as 



(e.g iLorimer fc Krameij |2005| ). where hem is the emission 
height and c the speed of light. 

The region bounded by the last open magnetic dipole 
field lines, and therefore possibly also the pulsar beam, is 
only circular if the magnetic axis is aligned with the rotation 
axis and for non-aligned rotators the polar cap is compressed 
in t he plain con taining the magnetic and the rotation axis 
(e.g. iBiggsll 19901 ). The beam shape becomes even more com- 
plex when the distortion of the mag netic field lines close to 
the li ght cylinder is considered (e.g. iRomani fc Yadigaroglul 
[l995|). Al so elongated beam shapes in t he latitudinal direc- 
tion (e.g. iNaravan fc Vivekanandll 19831 ) as well as an hour- 
glass shape for the beam of the binary pulsar B1913-f 16 
(jWeisberg fc Tayloill2002l ) have been proposed. 

Observationally, one measures the pulse width (A0) 
which is a function of p and the geometry of how the line 
of sight cuts through the beam . Und er the assumption that 
the beam is circular, iGil et al.l ([198J) derived 
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where a is the angle between the magnetic axis and the 
rotation axis and /3 is the impact parameter, the angle be- 
tween the line of s i ght a nd the magnetic axis at its closest 
approach. iRankinI ( 1990|) found that p for core emission is 
correlated with P~^'^ by using pulsars with interpulses for 
which a ~ 90°. Using this correlation in combination with 
the longitude dependence of the position ang le of the linear 
polarization (|Radhakrishnan fc Cookelll969l ) one can solve 
a and /3 for any pulsar with core emission. Several authors 
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used this method to derive a s imilar correla t ion between p 
and P for conal emi ss ion (e.g . Rankinlll993l : iKramer et al.l 
1 19941 : iGil et al.l[l99^ . iGouldl (|l994h found the same corre- 
lation by using the argument that the most narrow profiles 
are expected for pulsars with a ~ 90°. The implication of 
this correlation in conjunction with Eq. [T] is that the emis- 
sion height and the fraction of the polar ca p which is active 
is ind ependent of the pulsar period (e.g. iMitra fc RankinI 
|200^ 

iTauris fc Manchester! l| 19981 ) analysed published values 
for a of several hundred pulsars, showed they are incon- 
sistent with a random a distribution and conclu ded that a 
decreases on a time-scale of ^ 10^ yr. In contrast. iMcKinnonI 
l| 19931 ) concluded that the a distribution is random, provided 
that the beam was compressed in the latitudinal direction. 

In this paper, the issue of beam evolution is approached 
from a different angle. We do not rely on polarization data 
or a p correlation to determine a and /?. Rather, we will 
construct a beam model based on the observed fraction of 
pulsars with interpulses and their period distribution and 
the observed P — At/) plane is used directly. The interpulse 
population, the pulse width measurements and the model 
are described in sections 2, 3 and 4 respectively. The re- 
sults are presented in section 5 and their implications are 
discussed section 6. 



2 INTERPULSES 

The population of pulsars with interpulses is a powerful 
probe of pulsar beam properties. In the case where an in- 
terpulse is seen, with a separation from the main pulse of 
~ 180° in rotational phase, it implies that a must be ~ 90° 
and therefore gives a good count of how many orthogonal 
rotators there are in the population. We therefore searched 
the literature to obtain a complete census of the interpulse 
population. For this search, we excluded all millisecond pul- 
sars and pulsars in globular clusters, because the opening 
angle of the beams of milli-second pulsars are believed to 
have a different P depe ndence to that of normal pulsars 
(e.g. iKramer et al.ll 19981 ) and the values for the spin period 
derivative (P) of globular cluster pulsars are contaminated 
by their acceleration in the gravitational field of the cluster. 
We started with the catalogue of Taylor et al. (1993) which 
tabulated 14 interpulse pulsars out of a total of 527 objects. 
We then examined all pulse profiles in the major surveys 
conducted since that time. 

Table [T] lists the relevant papers, the number of pul- 
sars reported and the interpulses detected. Table [2] lists 
the population of pulsars with interpulses. The table in- 
cludes 3 weak interpulses not previously detected. Their 
profiles will be p resented in a forthcoming publication 
l|Weltevrede fc Joh nston 20Q3). In total there are 27 pul- 
sars with interpulses from a total sample of 1487 pulsars, 
i.e. 1.8% of the population. We note that the list should 
not be considered as complete. Some objects in the list may 
not be orthogonal rotators but are rather aligned rotators 
with wide b eams (see the debate for P SRs B0950-I-08 and 
B1929-I-10 in lEverett fc WeisberdboOll for example). On the 
other hand we may have missed some interpulse objects be- 
cause the flux in the interpulse is too weak to be detected. 

If the a distribution of pulsars is random, then of the 



Paper 



Ref. Np 



Nip 



Tavlor et al. (1993) 


(1) 


527 


14 


Foster et al. (1995) 




14 









19 





D'Amico et al. (1998) 


(2) 


84 


1 


Edwards et al. (2001) 




61 





Manchester et al. (2001) 




99 





Morris et al. (2002) 


(3) 
(4) 


119 
200 


3 


Kramer et al. (2003) 


2 


Hobbs et al. (2004) 


(5) 


178 


2 


Jacobv (2005) 




19 





Lorimer et al. (2005) 




8 
14 





Bureav et al. (2006) 





Lorimer et al. (2006) 


(6) 


131 


2 


Manchester et al. (2006) 




14 





Weltevrede & Johnston (2008) 


(7) 





3 



Total 



1487 



27 



Table 1. The used literature to find interpulses is listed in the 
first column. The reference codes are the same as in Table |2] 
Afpuiaars and A^ip are the total number of newly discovered normal 
pulsars in the discussed paper and the number of which have 
interpulses. 



JName 



BName 



Pis) 



Ref. 



J0534-t-2200 


B0531+21 


0.0331 


4.2 


X 


10' 


-13 


(1) 


J0826-I-2637 


B0823+26 


0.5307 


1.7 


X 


10- 


-15 


(1) 


J0828-3417 


B0826-34 


1.8489 


1.0 


X 


10- 


-15 


(1) 


J0834-4159 




0.1211 


4.4 


X 


10- 


-15 


(4) 


J0905-5127 




0.3463 


2.5 


X 


10- 


-14 


(7) 


J0908-4913 


B0906-49 


0.1068 


1.5 


X 


10- 


-14 


(1) 


J0953-f0755 


B0950-f08 


0.2531 


2.3 


X 


10- 


-16 


(1) 


J1057-5226 


B1055-52 


0.1971 


5.8 


X 


10- 


-15 


(1) 


J1126-6054 


B1124-60 


0.2027 


2.8 


X 


10- 


-16 


(7) 


J1302-6350 


B1259-63 


0.0478 


2.3 


X 


10- 


-15 


(1) 


J1549-4848 




0.2883 


1.4 


X 


10- 


-14 


(2) 


J1637-4553 


B1634-45 


0.1188 


3.2 


X 


10- 


-15 


(7) 


J1705-1906 


B1702-19 


0.2990 


4.1 


X 


10- 


-15 


(1) 


J1713-3844 




1.6001 


1.8 


X 


10- 


-13 


(4) 


J1722-3712 


B1719-37 


0.2362 


1.1 


X 


10- 


-14 


(1) 


J 1739-2903 


B1736-29 


0.3229 


7.9 


X 


10- 


-15 


(1) 


J1806-1920 




0.8798 


1.7 


X 


10- 


-17 


(3) 


J1808-1726 




0.2410 


1.2 


X 


10- 


-17 


(6) 


J1825-0935 


B1822-09 


0.7690 


5.2 


X 


10- 


-14 


(1) 


J1828-1101 




0.0721 


1.5 


X 


10- 


-14 


(3) 


J1843-0702 




0.1916 


2.1 


X 


10- 


-15 


(5) 


J1849-I-0409 




0.7612 


2.2 


X 


10- 


-14 


(6) 


J1851-I-0418 


B1848+04 


0.2847 


1.1 


X 


10- 


-15 


(1) 


J1852-0118 




0.4515 


1.8 


X 


10- 


-15 


(5) 


J1913+0832 




0.1344 


4.6 


X 


10- 


-15 


(3) 


J1932-I-1059 


B1929+10 


0.2265 


1.2 


X 


10- 


-15 


(1) 


J1946-fl805 


B1944+17 


0.4406 


2.4 


X 


10- 


-17 


(1) 



Table 2. The 27 interpulses of normal pulsars which can be found 
in the literature. The reference codes are the same as in Table [T] 



order of 5% of the pulsars ca n be expected to have an in- 
terpulse (e.g. iGil fc Hanlll996h . The measured percentage is 
lower, indicating a potential problem for such a model. Not 
only is the percentage of pulsars with interpulses an im- 
portant clue for beam models, more important still is their 
period distribution. This distribution is shown in Fig. [I] as 
a cumulative distribution (thick solid histogram). The dis- 
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Period 

Figure 1. The histogram (thick solid line) shows the observed 
cumulative distribution of pulse periods of the normal (i.e. non- 
millisecond) pulsars with interpulses. The lowest (thick) solid 
curve shows the cumulative period distribution of all normal pul- 
sars. The dotted curve shows the predicted distribution for the 
pulsars with interpulses in the case of a random a distribution 
(both for circular and non-circular beams). The dashed and solid 
curve are for a model including alignment of the magnetic axis 
for non-circular and circular beams respectively. 



tribution indicates the fraction of pulsars with interpulses 
y which have a period less than x. Compared with the pe- 
riod distribution of all normal (i.e. non-millisecond) pulsars 
(lowest thick curve of Fig. [l]) there are relatively many more 
short period pulsars with interpulses. In this paper we will 
try to construct a geometrical model which can account for 
the observed statistical properties of the population of pul- 
sars with interpulses. 



3 THE OBSERVED PULSE WIDTHS 
DISTRIBUTION 

We have embarked on a long-term timing campaign at the 
Parkes telescope in Australia to monitor a large sample of 
young pulsars with high spin-down energy loss rates E. This 
sample of pulsars is compared with archival Parkes data, 
which is published in various papers, in order to determine 
the differences between pulsars with high and low values for 
E (Wcltevrcde fc Johnston 2008 ). We refer to that paper for 
the details of the observations and data reduction. 

The dataset is used to determine the correlation be- 
tween the pulse width A(j) and P. Fig. [5] suggests that the 
slope in the Acp — P plane is less steep than —1/2. The 
slope measured by minimizing the x^ in log-log space (as- 
suming equal weights) is —0.31 ±0.05. As pointed out in the 
introduction, p is expected, and found, to be proportional 
to P~^'^ . An identical correlation can be expected in the 
A(j) — P plane, provided that a is independent of P. This 
measurement is therefore of interest to us because it can be 
interpreted as evidence for a P dependence of a. 

One could argue that the scatter of the points around 
the correlation in Fig. [2] is so large that a power law with 
a slope of —1/2 would fit the data equally well. We there- 
fore also measur ed the slope using the measured 10% widths 
at 1400 MHz bv lGould fc Lvnd i|l998D . excluding the milli- 
second pulsars, as a consistency check. Although there is 
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Figure 2. The measured profile 10% widths versus P at an 
observing frequency around 1390 MHz llWeltevrede &: JohnstonI 
200q). Pulsars with a S/N < 30 or those with substantial scatter 
broadening were excluded leaving 205 measurements. The solid 
line is a power law fit through the data, which has a slope of —0.31, 
comparable with the fit obtained from the data of lGould fc Lvnd 

ham)- 



some overlap in our samples, that of lGould fc Lvnd (|l998f ) 
contain many northern pulsars which are not pre sent in our 
sample. The slope measured in the iGould fc L vnc (19981) 
data is also flatter then expected (— 0.23±0.06), which makes 
us more confident that the observed trend is real. 



4 THE MODEL 

We attempt to explain three observational conclusions by 
creating a model which is as simple as possible. If one un- 
derstands the geometry one should first of all be able to 
construct a model which can predict the observed fraction of 
pulsars with an interpulse (1.8% or lower, see section 2). Sec- 
ondly, the model should be able to reproduce the observed 
period distribution of pulsars with interpulses. Thirdly, the 
model should explain why the observed A(/!> — P correlation 
is fiatter than the expected P~^" correlation. The effects 
of a possible alignment timescale of the magnetic axis and 
non-circular beams will be explored. 

The basic approach of this model is to synthesize a pop- 
ulation of pulsars in which the P — P distribution is identical 
to the observed distribution. We then assign a random ge- 
ometry to each pulsar and the subpopulation of pulsars with 
interpulses is identified using different model assumptions. 
The statistical properties of the synthesized population of 
pulsars with interpulses is then compared with the observed 
population in order to determine which model is able to re- 
produce the data best. 



4.1 The period distribution 

The model should first of all reproduce the observed P — P 
distribution for the total observed population of pulsars 
(those with and without interpulses), which is achieved by 
simply using the observed P — P combinations. Not only 
does this ensure that the resulting distributions are realis- 
tic, it also avoids the necessity to make an full evolutionary 
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model which would require additional assumptions about 
pulsar birth properties and their evolution. Note that we 
are therefore not making any assumptions about the parent 
distribution that gives rise to the observed distribution of 
P — P combinations. In the next subsection we will describe 
geometrical factors which could make the period distribu- 
tion of the pulsars with interpulses different to the total 
population of observed pulsars. 

The catalogue of p ulsars maintained by the ATNFu 
iJManchester et al.ll2005l ) was used to obtain the observed 
P — P values. Because the interpulses of milli-second pul- 
sars are not considered in this paper, all pulsars with an 
estimated surface magnetic field strength below 10^^ Gauss 
or a period faster than the Crab were excluded. PSR J 1808- 
1726, with a magnetic field strength of 5.4 x 10^" Gauss, is 
excluded from the interpulse statistics for consistency. The 
globular cluster pulsars have unreliable P values and are 
therefore also not considered. Finally, all the pulsars with 
P > 4 seconds were removed to exclude the catalogued soft 
gamma-ray repeaters and the anomalous X-ray pulsars. 



4.2 The beam model 

The first step is to draw random P ~ P pairs from the to- 
tal population of observed pulsars as described in Sect 4.1. 
Secondly, random values are generated for ct and 13. The 
magnetic pole and the line of sight are taken to be a ran- 
dom points on a sphere. This implies that the distribution of 
Q and C, = ot + j3 are sinusoida L The magn etic axis is allowed 
to align on a timescale Taiign (|Joneslll97q l 



a{t) = ao exp (-t/raiign) , 



(3) 



where ao is the random a value at birth and t is the age of 
the pulsar, which is taken to be the characteristic age. The 
final step is to determine if one or both beams of the pulsar 
intersect the line of sight given by the randomly generated 
geometry, which depends on the assumed beam shape. As 
pointed out in the introduction, different shapes of the polar 
cap have been proposed. Here we consider only the shape of 
a polar cap which is bounded by the last open magnetic 
dipole field lines in order to get a feeling for the effect of 
non-circular beam shapes. Such a beam is compressed in 
the plain containing the magnetic and the rotation axis and 
is r oughly elliptic. T he axial ratio S{q) of the ellipse is given 
by (|McKinnonlli993 ) 



£{a) = cos 5-y/cos [a ~ 5), 

where 5{a) follows from immerically solving 

2 tan ((5) = tan {a ~ 5) . 



(4) 



(5) 



The axial ratio £{a) varies from 1 (for an aligned beam) 
to ~ 0.62 (for an orthogonal beam). Circular beams can be 
assumed in the model by forcing £ = 1 for all ct. 

Apart from the beam shape, a beam size has also to be 
assumed. For the hal f opening ang le of the beam we use the 
relation measured bv lGouldl (|l994l ) 



p = 5.4°P" 



-1/2 



(6) 



which is consisten t with the re l ation reported by various 
other authors (e.g. lRankin|[l990l . 1 19931 : iKramer et al.lll994l : 
ICil et al.ll 19931 1. For non-circular beams, p of Eq. [S]is taken 
to be the half opening angle in the longitudinal direction. 
The conditions for the two pulsar beams to intersect the line 
of sight are 



\l3\^£{a)p 
and 
|180° -2«-/3| < f(«)P- 



(7) 



(8) 



http://www.atnf.csiro.au/research/pulsar/psrcat 



Random values of a and 13 are generated until at least one 
of these conditions are met. The width of the pulse profile 
can then be calculated with Eq. [2]for circular beams. 

As noted in section 2, the list of interpulses is not com- 
plete for two reasons. The first reason is that some inter- 
pulses might have been missed because the flux in the inter- 
pulse is too weak to be detected. Extremely weak interpulses 
could indicate that the pulsar beam just grazes the line of 
sight. Because the observed p (Eq.[S| has been derived con- 
sidering the 10% pulse widths, beams which just graze the 
line of sight are not counted as detections in the model (Eqs. 
[7] and O . This means that extremely weak interpulses are 
not accounted for in the model, as well as possibly in the 
observations. However, some of the observed interpulses will 
be wide beams of aligned rotators. Therefore the observed 
interpulse fraction should be considered as an upper limit. 



5 RESULTS 

Let us start with considering the simplest model, i.e. a ran- 
dom a distribution without alignment (raiign = oo) and cir- 
cular beams (£"(«) = 1). In this scenario our model predicts 
that 4.4% of the pulsars should have an interpulse, a fac- 
tor of at least two higher than the observed fraction. This 
model also clearly fails to explain the observed period dis- 
tribution of the pulsars with interpulses as the dotted curve 
of Fig. [T] lies below the observed distribution. This means 
that there are more short period pulsars with interpulses 
then predicted. Short period pulsars are expected to have a 
larger probability to show an interpulse because their beams 
are wider (Eq.|6]), which is the reason why the predicted dis- 
tribution of pulsars with interpulses (dotted line) lies above 
the distribution of all pulsars (lowest thick curve) . Neverthe- 
less this effect is not enough to explain the observed period 
distribution of the pulsars with interpulses. The predicted 
pulse width for this model is correlated with p-^-^^ ^ very 
close to the what is expected from Eq. (6] This model there- 
fore fails to explain the three observational properties out- 
lined in section 4. 

In order to make the model more accurately reflect the 
data it is necessary to produce fewer pulsars with inter- 
pulses, which can be done by relaxing the condition that 
forces the beams to be circular. Because the minor axis of 
non-circular beams is in the plane containing the magnetic 
and rotation axis, it is less likely that the beams of both 
magnetic poles intersect the line of sight. This in contrast 
to elongated beams, which will increase the fraction of pul- 
sars with interpulses and are therefore not considered further 
here. Indeed, when applying Eq. |4j the predicted fraction of 
interpulses is decreased to 2.3%. This is still a bit too high. 
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especially as the observed fraction is an upper limit. Because 
the ellipticity of the beam is independent of P it follows that 
the predicted pulse period distribution of pulsars with inter- 
pulses and the pulse width distribution are identical to the 
model with circular beams. Therefore the model with non- 
circular beams also fails to explain the observations. 

We next consider the effect of alignment of the pulsar 
beam. By setting Taiign = 7 x 10^ years (assuming circu- 
lar beams) the predicted fraction of pulsars with interpulses 
drops to the observed value. Also the predicted period dis- 
tribution of the pulsars with interpulses (thin solid line Fig. 
[l]) now agrees well with the observations, indicating that the 
model successfully predicts the observed under-abundance of 
long period pulsars with interpulses. This can be understood 
because long period pulsars tend to have larger character- 
istic ages, which means that they tend to be more aligned, 
which reduces the probability that the beams of both mag- 
netic poles intersect the line of sight. Finally, because the a 
distribution depends on P in this case, the predicted corre- 
lation between A(/!> and P will deviate from —1/2. The slope 
of the correlation becomes flatter (the slope is —0.40), as is 
observed. 

Finally, one can consider a model which both includes 
the ellipticity of the beam and alignment. In that case the 
timescale of alignment is much longer (raiign = 2 x 10^), 
because both effects reduce the predicted number of inter- 
pulses. As a consequence, the predicted period distribution 
of the pulsars with interpulses (dashed line Fig. [T| does not 
deviate much from the model without alignment and there- 
fore does not fit the data as well as the model with circular 
beams. 

If pulsar beams align over time, it would imply that the 
total population of pulsars has an a distribution which is 
skewed to low values compared with a random distribution. 
This affects the average beaming fraction because i t is a 
function of a (e.g. Eq. 7 of i Tauris fc Mancheste3ll998l ). The 
beaming fraction is the fraction of the celestial sphere swept 
out by the beams of a pulsar, hence it is the probability that 
a pulsar is observable for a random line of sight. Aligned 
beams are less likely to intersect the line of sight and are 
therefore less likely to be observed. The average beaming 
fraction for the observed population of pulsars, assuming 
circular beams, is 8.4% and 17% for the model with and 
without a evolution respectively. 

For a random a distribution the fraction of pulsars with 
beams which continuously intersect the line of sight is tiny 
(0.01%), because of the a dependence of the beaming frac- 
tion in combination with the effect that aligned pulsars are 
less likely to be formed in the first place. This means that ba- 
sically no pulsars with extremely wide profiles are expected 
in the known population of pulsars. This fraction appears 
too low, as we have already argued that some of the inter- 
pulses in Table [5] originate from a single magnetic pole. A 
more realistic fraction is predicted by the model assuming 
circular beams including a evolution (1.6%), while Taiign is 
too long to make a difference in the case of non-circular 
beams. This can be seen as additional evidence for circular 
beams and a significant effect of a evolution. 



6 DISCUSSION 

We have shown that the observed fraction of pulsars with 
interpulses, their period distribution, the fraction of pulsars 
with extremely wide profiles and the observed pulse width 
versus pulse period correlat ion (as measured in our data and 
that of lGould fc Lvnelll998l ) is inconsistent with a random a 
distribution, even when non-circular beams are considered. 
The main problem is that too many pulsars with interpulses 
are predicted, especially those with long periods. A number 
of explanations are possible. 

First, one could assume that long period pulsars are less 
likely to have two active poles. Not only is this assumption 
ad-hoc, it also does not explain the observed correlation be- 
tween the pulse width and the pulse period. Our model is 
kept simple and intuitive by assuming that both magnetic 
poles always emit radio emission, the beams are directed in 
exactly opposite directions and have equal widths and lumi- 
nosities. We therefore aim to find a geometrical solution for 
the problem. 

Secondly, one could assume a different powerlaw rela- 
tion between p and P. In order to correctly predict the num- 
ber of pulsars with interpulses and their period distribution 
the observed powerlaw relation (Eq. [6]) should be steeper 
and the constant smaller. This could physically mean that 
the fraction of the polar cap which is active or the emission 
height is P dependent. Not only is such a model incompati- 
ble with the observations from which Eq.|6]is derived, it also 
results in a highly unrealistic pulse width distribution (the 
predicted pulse widths are too narrow and the correlation 
with P is too steep). 

Thirdly, the pulsar beam may not be circular, which 
is expected if the polar cap is bounded by the last open 
magnetic field lines. Although elliptical beams which are 
compressed in the plain containing the magnetic and the 
rotation axis could help to reduce the number of predicted 
interpulses, it is not enough. Moreover, the ellipticity is not 
expected to be P dependent. This implies that this model, 
like that for circular beams, does not correctly predict the 
period distribution of the pulsars interpulses and the pulse 
width period correlation. The fraction of pulsars with ex- 
temely wide profiles is also too low. 

The only geometrical effect left, as far as we can see, to 
explain the observations is to relax the assumption that a 
is random. In order to explain the observed fraction of pul- 
sars with interpulses the a distribution should be skewed to 
low values. The distribution for short period pulsars should 
be skewed less than that of long period pulsars in order to 
explain the observed pulse period distribution of the pul- 
sars with interpulses. The skewness therefore appears to be 
caused by an evolution in time of a rather than by a non- 
random birth distribution which is fixed in time. 

In the previous section it has been shown that the model 
assuming circular beams can explain the observations well 
when the magnetic axis is allowed to align from a random 
distribution at birth with a timescale of '^ 7 x 10^ years. The 
fact that the fraction of pulsars with interpulses, their pe- 
riod dependence, the fraction of pulsars with extremely wide 
profiles and the pulse width distribution can be explained 
by adding just one extra model parameter is encouraging. 
Moreover, by allowing a to evolve it is not necessary to aban- 
don the observed p — P relation (Eq.[6]) or to resort in ad- hoc 
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assumptions about magnetic poles which are not always ac- 
tive. 

The timescale is comparible with the timescalc derived 
from the observed a distribution by iTauris fc Manchester 
ljl99a ). Measuring a is far from trivial and for most pulsars 
it is, at best, only poorly constrained. It should be stressed 
that although we discuss a evolution, we do not rely on 
(indirect) measurements of a. Therefore our results are an 
important addition to the existing evidence for alignment. 

The pulse period distribution of the pulsars with in- 
terpulses can be used to distinguish between beam shapes. 
The data show that the derived distribution assuming cir- 
cular beams fits the data better than that assuming ellip- 
tical beams. This suggests that the beams are more cir- 
cular than expected from a simple model which has the 
polar cap bounded by the last open dipole magnetic field 
lines. For example, more complicated models which take 
into account the distortion of the magnetic field lines close 
to the light cylinder predict roughly circular beams (e.g. 
iRomani fc Yadigaroglulll995l ). 

Our estimation of Taiign should be considered to be 
very rough because of a number of reasons. First of all, the 
timescale is based on the observed characteristic age distri- 
bution of pulsars. The characteristic age should be consid- 
ered to be only a rough estimator of the true age, because 
the details of pulsar spin-down are not well understood. Sec- 
ondly, the a evolution could have a different form than the 
assumed exponential decay. For instance one could imagine 
that the decay rate is a function of a. Thirdly, Taiign also 
depends on the shape of the radio beam. If the beam is non- 
circular then Taiign should be larger in order to explain the 
observations. Finally, the timescale is based on the observed 
fraction of interpulses. It is assumed that all the observed 
interpulses are caused by emission from the opposite pole, 
which is unlikely to be the case. Therefore the observed frac- 
tion of pulsars with interpulses should be seen as an upper 
limit, resulting in an overestimation of Taiign. 

Although our estimated value of Taiign is very rough, the 
observations can only be explained if the effect of alignment 
is important and therefore happens on a timescale compara- 
ble to the age of a typical pulsar. Beam evolution has some 
important consequences for population studies, evolutionary 
models and pulsar searches, for two reasons. 

The first is that the beaming fraction depends on Taiign 
such that aligned beams are less likely to intersect the line of 
sight. Therefore a consequence of alignment of pulsar beams 
over time is that older pulsars become harder to find, hence 
the inferred total population of old pulsars should be larger. 
The average beaming fraction for the observed population 
of pulsars, assuming circular beams, is 8.4% and 17% for the 
model with and without a evolution respectively. Therefore 
the inferred total population of pulsars is about twice as 
high when a evolution is considered. 

Secondly, one can ask the more fundamental question 
of what alignment implies for the spin-down of pulsars. The 
evolution of a must be the result of the braking torque, 
hence it is related to the spin-down of pulsars. If the torque 
changes a, it seems reasonable that the torque also depends 
on Q. It is therefore not unlikely that the pulsar spin-down 
itself has an a dependence. The spin-down of a rotating 
dipole field in a vacuum depends on a (_B oc sin^ a, e.g. 
|jacksonlll962l ). but how different this is in the presence of a 



pulsar magnetosphere is unclear. An a dependent spin-down 
implies that the magnetic field strength, the characteristic 
age and spin-down energy loss rate of a pulsar also depends 
on a, which has important consequences for evolutionary 
models and population studies. 

We have shown that the observed population of pulsars 
with interpulses places strong restrictions on beam models. 
Future instrumentation, such as the Square Kilometre Ar- 
ray, will greatly enlarge the observed population of pulsars, 
enabling the possibility to distinguish between models with 
different beam shapes and different functional forms for the 
Q evolution. This will contribute to a better understanding 
of how breaking torques in neutron stars work. 
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